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We study electronic transport through an Aharonov-Bohm interferometer with single-level quan- 
tum dots embedded in the two arms. The full counting statistics in the shot-noise regime is calculated 
to first order in the tunnel-coupling strength. The interplay of interference and charging energy in 
the dots leads to a dynamical channel blockade that is tunable by the magnetic flux penetrating 
the Aharonov-Bohm ring. We find super-Poissonian behavior with diverging second and higher 
cumulants when the Aharonov-Bohm flux approaches an integer multiple of the flux quantum. 
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I. INTRODUCTION 

The study of full counting statistics (FCS) of charge 
transport through mesoscopic systems has become a 
well established field. A number of theoretical ap- 
proaches for calculating the cumulant generating func- 
tion have been developedfii^i^i^i^i^ extended to inter- 
acting systems,— i^iSiii^ and memory effects have been 
included F^iii^iii^ Experimentally, it has become possible 
to count single electrons in real time as they pass through 
a system of quantum dots i^^'^^i^^ Despite the detector 
performing a projective measurement, interference has 
been observed in specially designed systemsji^, Apart 
from the conccptionally straightforward idea of counting 
individual electrons there have been numerous propos- 
als for counting statistics detectors, involving qubits and 
superconducting systemsJ^'i^'^SiiSiiS^ 

A useful reference to compare transport statistics with 
are Poissonian processes, which describe uncorrelated 
events of charge transfer. Correlations between the trans- 
port events change the transport statistics. For fermionic 
systems, correlations typically lead to a reduction of the 
current noise. However, there are also various scenar- 
ios in which the noise is increased. A prominent example 
that exhibits super-Poissonian noise is electron bunching, 
in which periods of high and low (or zero) current alter- 
nate. Bunching behavior has been discussed in complex 
geometries such as beam splitters^^ and serial doubleSiiS 
and triple^^ quantum dots. It may be understood as a 
consequence of a system's bistabilityj ^^i^^ A further ex- 
ample for a bistable system is a quantum shuttle, which 
exhibits enhanced noise at the transition from the tun- 
neling to the shuttling regime 

Super-Poissonian transport behavior can already be 
found in a single quantum dot, e.g. when lifting of 
the spin degeneracy of the level results in different tun- 
neling rates for the two spin statesi^ A similar effect 
was measured in a quantum dot with two states cou- 
pling differently to the leads either due to their differ- 
ing spatial extensioni^ or due to spin-dependent tunnel- 
ing to ferromagnetic leadsi '^^'^^ Transport through a sin- 
gle level can also exhibit bunching if a second level in 
its vicinity interrupts transport by means of Coulomb 
interactioni2^i2ii^i2& In a similar way, driven transitions 



between two levels may give rise to enhanced noise 

Also the counting statistics of transport through par- 
allel double quantum dots may exhibit super-Poissonian 
behavior. Enhanced noise of the cotunneling current has 
been found in the presence of ferromagnetic leadsi^ In 
DQDs with normal leads enhanced noise was predicted 
for spinless electrons witb^^i^Si^ and without^ inter-dot 
Coulomb interaction and for spinful electrons in dou- 
ble dots with an inter-dot tunnel coupling that leads to 
a splitting between symmetric and antisymmetric one- 
electron states in the double dot4i 

In this paper we study a parallel double quantum dot 
(DQD) as shown in Fig. [1] The dots are weakly coupled 
to two leads, operated in the shot-noise regime, with an 
Aharonov-Bohm flux <!> enclosed by the two paths. We 
assume the two dots to be sufficiently separated such that 
there is no direct tunneling and no inter-dot charging 
energy. On the other hand, we assume a strong Coulomb 
interaction within each dot. The system can, therefore, 
accommodate at most two electrons. The spin degree 
of freedom will turn out to be a vital ingredient, since it 
necessitates the description of the doubly-occupied states 
by singlet and triplets. These will turn out to be critical 
for the decomposition of the system's Liouville space into 
two disjoint parts, resulting in diverging cumulants. 

The possibility to entangle the spins in the dots is used 
in the context of quantum information processing4^i^ 




FIG. 1: The double-dot Aharonov-Bohm interferometer con- 
sists of two quantum dots {u and d) connected to two leads in 
parallel. The paths through the dots enclose a magnetic flux 
"I>, giving rise to Aharonov-Bohm interference. 
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It was recently pointed ou t ' that in non-equilibrium 
situations an imbalance between spin singlet and triplet 
states in the DQD can be generated, based on a scheme 
that is similar to coherent population trapping4ii^ In 
contrast to double dots with direct tunneling between 
the dots, the symmetric and antisymmetric one-electron 
states remain energetically degenerate for weak tunnel 
coupling to the leads. The same holds true for the singlet 
and triplet two-electron states. The imbalance between 
singlet and triplet is, therefore, a genuine non-equilibrium 
effect that relies on the interplay between coherent tunnel 
coupling and Coulomb interaction4^i^ 

In the Aharonov-Bohm interferometer discussed in 
this paper, the singlet-triplet imbalance depends on the 
Aharonov-Bohm phase, which can be tuned by the mag- 
netic flux enclosed by the interferometer arms. Its origin 
is related to the fact that the source and drain electrons 
only couple to certain linear combinations of the upper 
and lower dot levels. Most dramatic effects are expected 
for Aharonov-Bohm fluxes that are close to integer multi- 
ples of the flux quantum. In this case, as we will discuss in 
more detail below, the singlet and triplet states belong to 
two different subspaces of the double dot's Hilbert space. 
These subspaces are nearly disconnected from each other 
and are described by different transport characteristics. 
As a consequence, we will find bunching behavior that 
leads to not only an enhancement but even a divergence 
of the second and higher cumulants as a function of the 
Aharonov-Bohm phase. 

The paper is structured as follows: In Section we 
specify the system model. Section IIIII describes how to 
calculate the full counting statistics. The results are dis- 
cussed in Sec. IIVI where we illustrate the mechanism of 
super-Poissonian statistics and discuss several limiting 
cases. In Sec. |V]we demonstrate the origin of the effect 
by comparing to a spinless model. Finally, we conclude 
in Sec. ED 



II. SYSTEM 

The double-dot interferometer shown in Fig. [T] is de- 
scribed by the Hamiltonian 

H = Hu + Hd + HL + HR + HT. (1) 

The quantum dots. Hi = J^a^i'^i'^a-i + Un-^iUii for 
i ^ u,d, are described as Anderson impurities with spin- 
degenerate electronic levels and charging energy U for 
double occupation. Throughout this paper we are inter- 
ested in the regime of strong Coulomb interaction {U 
greater than all other energies), so that at most sin- 
gle occupation of each dot is allowed. Furthermore, 
we concentrate on the situation when both quantum- 
dot levels are tuned close to each other. We define 
the average level energy as e = (e^ + ed)/2. Each of 
the leads is described as a reservoir of non-interacting 
fermions Hr = J^kaa^rka '^Ika^rka with iudiccs for lead 



r G {L,i?}, momentum k and spin a. The tunneling 
Hamiltonian Ht = i HT,ri consists of parts for tun- 
neling between each dot i and each lead r, 

^T." = t„alka^„, + H.C., (2) 

k.a 

with the ffux dependence included in the phases of the 
tunnefing amplitudes tL,d = liije"^/^, tL,u = liije""^/'', 
tR..u = l^jjle"^/", and tR^d = l^flle"*"*/", according to 
Fig. [T] The phase (f> is related to the magnetic fiux <& 
through the ring as = 27r<I>/<I>o, where = h/e is the 
ffux quantum. The tunneling rate through interface r 
is quantified by Tr/fi ~ 2tt \tr\^ / h. For simplicity, we 
assume the density of states pr and the tunneling am- 
plitudes tr to be independent of energy, which implies 
constant tunneling rates. 

In addition to strong on-site Coulomb repulsion, 
we assume no interdot interaction, so that the entire 
system can be occupied by at most two electrons. 
The probabilities to find the system empty and singly 
occupied are pq and pi, respectively. Charging the 
empty system with an electron of spin a from the left 
lead results in the state (e^^ct^ -I- e~*$ c^^)/-\/2 |0) = 

(e*$|cr,0) + e-'i \0,a))/V2 = \+) ^. This state is 
not fully described by the probability of single occu- 
pation pi ^ Pu + Pd, it rather needs to be further 
specified by off-diagonal elements of the density matrix 
Pv ~ (W) il^l) where fi and ly label the nine dot states 

|o),|t,o),U,o),|o,T),|o,i),|t,T>,li,i),IT,i>,li,T>. 

We summarize them as an isospin in the two- 
dimensional Hilbert-space of the two dot levels 

la = {Ia,x,Ia,y,Ia,z) = {p'i + pt iip"^ - pfj , Pu " Pd) /2. In 

this basis the state reached by tunneling in from the left 
lead is isospin-polarized along = (cos 0/2, sin (/)/2, 0). 
The right lead is correspondingly isospin-polarized along 
rii^ = (cos 0/2, — sin 0/2, 0). 

Due to strong Coulomb interaction on the dots double 
occupation of the system is allowed only if one electron is 
found in each dot. This means that by sequential filling 
from the source lead only the singlet state \S) = (IT^ i) ^ 
|i,t))/V2 is accessible. The three triplets |r+) = |t,t), 
|T_) = U, i), |To) = (It, I) + li, T))/V2 can be accessed 
only indirectly as we will see later. Spin symmetry of 
the Hamiltonian requires all triplets to be occupied with 
equal probability pT+ — Pt_ = PTo ~ Pt in the stationary 
limit. While in principle the density matrix may contain 
9x9 = 81 elements, this number reduces to the following 
seven: po,pi,ps,PT and three isospin components I (spin 
symmetry requires the occupation of all triplets to be 
equal). 



III. MASTER EQUATION AND PCS 

The time evolution of the system's density matrix can 
be described by a number-resolved generalized master 
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equation 
d 



dt 



(3) 



= E dt'j2w^^i::{N-N\t-t')p^;iN',t'), 

where p'^{N,t) is the component of the density matrix 
under the condition that N electrons have passed the 
system after time t. In general, the energy difference 
{e^ — Eu) appears in the presence of off-diagonal matrix 
elements. It, however, vanishes for symmetric dot level 
energies, e„ = e^. 

The kernel of this equation can be obtained using 
a diagrammatic real-time technique formulated on the 
Kcldysh contour. It allows for a perturbative expansion 
in the coupling strength, which we abort in lowest order 
r to describe the weak-tunneling limit. For a detailed 
derivation of this diagrammatic language and its rules for 
the calculation of diagrams we refer to Refs. i45,46,49,53. 
Not described in these references is the inclusion of the 
counting field x^ the Fourier-conjugated variable of the 
transferred charge N . It is introduced at each junction 
by replacing the tunnel matrix elements in the Hamilto- 
nian as tr tre^^^'^l'^ with xl = — Xi? = x/2, where 
the positive (negative) sign is taken for vertices on the 
upper (lower) branch of the Kcldysh contour. The count- 
ing field is thus attached to both interfaces in such a way 
that only electrons passing the entire system contribute 
to the statistics. 

The cumulant generating function is defined as 



S{xM)) = In 



.N=-oo 



(4) 



where P{N,tQ) = I]^P^(A^,to) labels the probability 
that N electrons have passed the system and are the 
diagonal elements of the systems density matrix, i.e. the 
occupation probability of state The cumulants can 
be obtained by taking derivatives with respect to the 
counting field k„ = (-i)"(e"7io)9"/9x"S'(x)- The den- 
sity matrix elements p'j[ are the solutions of the master 
equation Eq. ([3]) in the steady state. In order to solve 
Eq. ([3]), we first Fourier-transform the equation with re- 
spect to iV, thereby introducing the counting field x and 
second perform a Laplace transform in time. We obtain 



the cumulant generating function by following the same 
steps as in Ref. [12, in which, however, only generalized 
master equations for a diagonal density matrix were con- 
sidered. This is, e.g., sufficient to describe Aharonov- 
Bohm interferometers that contain one quantum doti^ 
In order to describe a double-dot Aharonov-Bohm inter- 
ferometer, however, we need to extend the approach of 
Ref. [l^ to include also non-diagonal density matrix ele- 
ments. 

It turns out that this extension is quite straight- 
forward. All the formal steps of Ref. [l2 remain the 
same. The only difference is that when writing the 
generalized master equation in a matrix notation, off- 
diagonal matrix elements need to be taken into ac- 
count as well, i.e., 7r(x) = W(x) • f^{x) where tt = 
(p^i , p^^ , . . . , pi'^"^ , , . . . , p^Z ) first collects all di- 
agonal and then all off-diagonal matrix elements of the 
density matrix (the indices fii,pi,ai label the system 
states.). As a consequence, we need to redefine the vec- 
tor = (1, . . . , 1, 0, . . . , 0), which is needed to com- 
pute P{N) = • 7r(iV). The central result of Ref. [H 
is still valid, namely Eq. (5) and its generalization to 
non-Mar kovian orders. It relates the cumulant generat- 
ing function to the eigenvalue A(x) of the kernel W(x), 
whose real part has the smallest absolute (negative) 
value. The only modification is that the matrix W(x) 
is enlarged since it allows for transitions from and to off- 
diagonal states as well. In lowest order in the tunnel 
coupling strength, the cumulant generating function is 
simply given by 



Six,to) = toXix)- 



(5) 



In many cases the matrix W is too large and complex 
to obtain its eigenvalues analytically. In such situations, 
one can nevertheless obtain all the cumulants recursively. 
This has been demonstrated for master equations similar 
to Eq. ([3]) by Flindt et al— This technique can be applied 
to systems with non-diagonal density matrix elements 
without modification. Details and analytic forms for the 
the first cumulants can be found in the appendix and the 
literature. 

We present the master equation for the case of degen- 
erate dot levels (i.e. level splitting Ae = 0) using an in- 
tuitive notation distinguishing between occupation prob- 
abilities p = {po,Pi,ps,Pt) and isospin I [in the above 
notation tt = (p,I)]. 



/ _4/^ e-*x-(l-/,) \ 



r=L,R 





V 



e^x- i/,, -2(1 - U) 
e*x-|/, -2(1-/.) j 



• p + r. 



/ e-^^'- 2(1 - /,) \ 
-2(1 - 2U) 

\ -e*^- 3/, j 



I Xlr (6) 



-I = ^ [e'^^ 2frpa + (1 - fr)pi + e-'x- (1 - fr)ps ~ e"''^- (1 - fr)pT] - r,(l + /,)I 



(7) 



r=L,R 
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In the following we will only be interested in the shot- 
noise regime eV ^ ksT. In particular, we consider the 
situation that the dot levels lie inside the energy window 
defined by the Fermi levels of the leads, so that /L(e) = 1, 

IV. FLUX DEPENDENT COUNTING 
STATISTICS 

To illustrate the origin of the singlct-triplct imbalance 
discussed in Ref. IH, we transform the basis of the double- 
dot states with a transformation matrix S such that 
Stt = (po, 2| -f I-ni,ps, ^-I-ni,pT,I - (n. x ni),4). 
The first three elements correspond to the double dot 
being empty |0), singly occupied in the symmetric state 
l+>L = {e"''^^\cr,0) + e-'"f'/'^\0,a))/V2, that is reached 
by tunneling in from the left lead, and doubly occupied 
with a spin singlet \S). In the following this set of states 



is referred to as the -I- subspace, named after the isospin 
in the singly occupied state. This basis choice is moti- 
vated by the fact that for flux cj) = 27rm, the + subspace 
is not connected by tunneling to the remaining states. 
These other states can again be divided into two uncou- 
pled subspaces. The first one consists of single occupa- 
tion in the state \-~)r = 

(this becomes the antisymmetric state for (f) = 27rm), 
and double occupation with a spin triplet |T). The set 
of these two states is referred to as the — subspace. The 
remaining components of the isospin that are orthogo- 
nal to rii are called the _L-subspace. In the case (p = 
these states are not populated and the system dynamics 
is governed entirely by the -I— and — subspaces. 

The basis change includes no approximation and con- 
tains the same information found in Eq. \6\ since states 
and rates were transformed together. In the new ba- 
sis the master equation assumes block-diagonal form for 
(t) = 0: 
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(8) 



where the following definitions were used: G± = r4(l± 

cos I), r = rL + r^j and r^o = 2rL + Tr. Note that 
for = the matrix assumes block-diagonal form, since 
G+ ^ and G- ^ 0. 

It turns out that for zero flux the master equation as- 
sumes block diagonal form and the +- and — subspaces 
decouple. The -|- subspace is found in the upper-left 3x3 
block. In presence of a flux it is coupled by tunneling, 
i.e. a change in the charge state, to the 2 x 2-dimensional 
— subspace located in the middle. This coupling cannot 
be described simply by a rate ri+^_. Instead there are 
six possible transition paths which occur with four differ- 
ent rates. The paths and rates can be read off from the 
2x3- and 3 x 2-blocks in Eq. ([8]) that contain factors G^ . 
The two rightmost columns and two lowest lines of Eq. [8] 
describe the intermediate occupation of the _L-subspace 
occurring for ^ 27rm. 

As for (j> = 27rm the system decomposes into two un- 
coupled subsystems it is no longer possible to calculate 
its counting statistics as described in Ref. [H. There are 
two independent stationary solutions of the master equa- 
tion ©-one for each subspace. This means that the (long 
time) counting statistics would unphysically depend on 
the initial condition. In realistic systems some kind of 



coupling would always be present, lifting the degener- 
acy. Thus we study only small but finite values of the 
flux (j) Si 27rm so that the full counting statistics is well 
defined. 



A. Channel Exclusion for (f> = 27rm 

The separation of the system's Hilbert space into two 
separate subspaces has consequences for the transport 
statistics: As discussed in the introduction, systems with 
several states differing in average current exhibit bunch- 
ing. The larger the current difference is, the more en- 
hanced is the noise. As the rate with which the system 
switches between the states is decreased, the noise is ex- 
pected to be enhanced further. In the system discussed, 
the coupling can be decreased to zero by tuning the mag- 
netic flux to 4> = 2TTm. In this case, the subspaces decou- 
ple and all normalized cumulants diverge as the effective 
charge goes to infinity (see Fig.[2|)i^ We emphasize that 
the divergence of the normalized cumulants is not caused 
by a vanishing current, whose dependence on flux is ap- 
proximately cosinelike (see Fig. [3]) . 

To specify the picture outlined above we analyze the 



5 







— n=2 

— n=j 
■ ■ ■ n=4 












— H h-T^^ 1 i-^-^H ^ — 

1 1 

1 1 1 1 1 



71/2 71 37t/2 27C 
AB-Phase (|) 



FIG. 2: (Color online) The normalized nth cumulants h„/ki 
diverge for 2nm in a symmetric system {Tl = T^j), 
due to the competition of two channels belonging to different 
system states. The width of the divergence is governed by the 
relaxation rate between these two states. 



The complex internal dynamics distinguish our system 
from others in the hterature in two ways: First, the two 
states are not just diflFering in current, but are each char- 
acterized by their own distribution function. Second, 
the transitions between the subsystems cannot just be 
described by simple rates r+„_. Instead there are six 
possible transition paths which occur with four different 
rates. These paths and rates can be found in the 2x3- 
and 3 x 2-submatrices in Eq. ([8]). 

Our work also differs from other studies of the counting 
statistics related to singlet and triplet statea^i^ in so 
far as these studied the different statistics of two-particle 
states propagating along a device, while in our case only 
a single particle propagates at a time. 



properties of the subspaces separately. For this purpose 
we calculate the cumulant generating functions associ- 
ated with the 3x3- and 2 x 2-subspaces for singlet and 
triplet, respectively. These are defined as the eigenvalues 
of the corresponding sub-matrices of Eq. ([8]) . 

The cumulant generating function for the — subspace, 
Eq. ^ is that of a two-state system. 



4(3rL)(2rfl) 

(3ri + 2Tny 



The transition rates are ST^ for filling and 2Tf( for emp- 
tying the double dot. This can be understood by count- 
ing the possible realizations of each state: If the dot is 
in a triplet state, taking away either of the two electrons 
results in single occupation. On the other hand, start- 
ing from a singly occupied state there are three triplets 
which can be accessed by tunneling into the system. 

The cumulant generating function for the + subspace 
is too complex to be shown here. It describes a three- 
state system, with transition rates as can be read off from 
Eq. ([51). Below, we will give compact analytic expressions 
in the limit of very asymmetric tunnel couplings to the 
left and right lead. 
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FIG. 3: The current through the double dot system is subject 
to Aharonov-Bohm oscillations. These are almost cosine-like 
so that the current remains finite for d) = 0. 



1. Asymmetric Tunnel Coupling 



In Fig. m we show the second, third, and fourth 
normalized cumulant for asymmetric tunnel coupling to 
source and drain, parametrized by the asymmetry pa- 
rameter a = {Tl — rfl)/r, where F = -I- F^j denotes 
the total coupling. For a ^ — 1, the bottleneck for trans- 
port is the tunnel barrier between source electrode and 
double dot. We find that the width of the divergence as a 
function of the Aharonov-Bohm flux is slightly increased. 
Furthermore, the divergence of the third cumulant has 
changed its sign as compared to the case of symmetric 
coupling. The cumulant generating functions for the +- 
and — subspace at = 27rm, 



-^F(a + l)V>^(e^>^-l), (10) 
= 3Fi(e^>^-l) 

^lT{a+lfe'^{e'^-l), (11) 



become Poissonian for a ^ — 1, with different tunneling 
rates for the -|— and — subspace. In addition, the proba- 
bility to find the double dot in a + subspace is higher than 
that for the — subspace, P+ =pQ+p+-\-ps = | 4- 0(0^) 
and P- ^ P-+PT = ■^+0{a^). These probabilities were 
obtained from the probabilities calculated in the non- 
degenerate case (p and then taking the limit 0^0. 

In the opposite limit, a 1, the width of the diver- 
gence is strongly suppressed. In this case, tunneling out 
of the double dot to the drain electrode defines the bot- 
tleneck of transport. At first glance a stronger bunch- 
ing may be expected, since with decreasing tunnel cou- 
pling to the drain the coupling between the -I— and — 
subspaces is reduced. However, it turns out that in the 
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FIG. 4: (Color online) The normalized nth cumulants 
of the non-symmetric system diverge (left plot; a — {Fl — 
^r)/{Tl + Vr) = -0.9, right plot: a = +0.9). The width 
of the divergence is slightly enhanced a — > — 1 but strongly 
reduced for a ^ +1. 




Asymmetry a=(rj^-rj^)/(rj^-Hrj^) 



FIG. 5: (Color online) For stronger coupling to the drain (a < 
0) the current associated with the + subspace is larger that 
that for the — subspace. Additionally, the occupations of the 
subspaces differ. This results in super-Poissonian statistics. 
For a — » 1 both currents and occupations are equal. This 
yields Poissonian statistics. 



limit a ^ +1. the cumulant generating functions 

_ir(a- l)3e2*x(e*'^ - 1) (12) 

-ir(a-l)Vx(e'>^-l) (13) 

become identical. Therefore, the transport statistics is 
the same for both subspaces, and bunching does not ap- 
pear anymore. Figure [5] shows the current /-|_ (/_) under 
the condition that the system is in the -I- (— ) subspace as 
a function of the the asymmetry. These currents are ob- 
tained from the cumulant generating functions for each 
subspace. They become equal for a —f +1. 



2. The influence of relaxation 

In experiments interaction with the environment can 
be expected to induce various relaxation mechanisms. 
Relaxation of the isospin, mediated by electric interac- 
tions, will be of particular importance. We model it by 
introducing a relaxation term with rate 17i in the master 
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FIG. 6: (Color online) Normalized cumulants in the presence 
of isotropic isospin relaxation (rate Q.i = F/IO). The diver- 
gent cumulants are suppressed: while the Fano-factor is sub- 
Poissonian for the relaxation rate shown, higher cumulants 
still show super-Poissonian behavior. 



equation Eq. ^ that reduces the isospin isotropically 
d 



dt 



(14) 



rol 



The effect of this relaxation is primarily a reduction of the 
visibility of the AB-signal, due to the electrons loosing 
their coherence. Furthermore, it leads to an effective 
coupling of the +- and — subspaces. Correspondingly the 
bunching effect is weakened, resulting in the cumulants 
assuming finite values also for = 2Tim (see Fig. (5]). The 
figure shows the situation for fij = F/lO. For sufficiently 
fast relaxation all cumulants become sub-Poissonian, but 
as can be seen from the figure, the higher the moment, 
the faster is the required relaxation rate. 

Another relaxation mechanism may be given by spin- 
flip processes converting singlets into triplets and vice 
versa. Since we summarized all the triplet occupations 
in pt, terms have to be added to the master equation in 
the following way: 



PS 



rol 



Pt 



-flSTPS + 3^stPt (15) 

-flsT PS - SrisT PT- (16) 



The factor of three is required to take into account that 
the triplet probability corresponds to three states, while 
there is only one singlet. 

The -I— and — subspaces are now directly coupled and 
the divergencies vanish more rapidly as a function of 
the relaxation rate than in the case of isospin relaxation 
(Fig. [7]). However, since such a relaxation is mediated 
magnetically, it can be expected to be much slower than 
isospin relaxation. 

A third way of coupling the -I— and — subspaces is not 
related to relaxation. As can be seen from the master 
equation, Eq. a detuning of the energy levels Ae = 

— Ed gives rise to additional terms. It turns out that 
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FIG. 7: (Color online) Normalized cumulants in the presence 
of singlet-triplet relaxation (rate ilsT ~ r/200). The influ- 
ence of the S-T relaxation is much stronger than that of the 
isospin relaxation. 
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FIG. 8: (Color online) Normalized cumulants in the presence 
of a finite level detuning Ae — l/SF. Symmetry of the flux- 
dependence about (j) = TT is lost, but enhancement of the 
cumulants are still clearly visible. 



for small detuning, Ae ~ F, this results in precession of 
the isospin about the axis n = (0, 0, 1), 



dt 



= Ae n X I. 



(17) 



Since three spatial directions (n^, and n) appear 
in the master equation, the symmetry of the flux- 
dependence about = TT is lost. In other words, the 
statistics depends on the direction of transport even in 
the case of symmetric coupling F^ — Tn. A similar effect 
was predicted for a quantum dot with three ferromag- 
netic leads.— For reasonably large values of the detuning, 
Ae = F/3, enhancement of the moments clearly persists 
(Fig.®. 



B. Super-Poissonian statistics for (j) = (2n -I- l)7r 

Figure|3]reveals that the statistics for (j> = (277-|-l)7r arc 
also peculiar for a — > — 1: The noise and the third nor- 
malized cumulant arc enhanced beyond the Poissonian 




FIG. 9: Internal system dynamics for flux (j) — tt and strong 
coupling to the drain a — > — 1. The double lines denote tran- 
sitions much faster (with rate 2F_r) than the single lines (with 
rates oc Fl). 



value, while the fourth normalized cumulant remains neg- 
ative. The enhanced noise can be understood by studying 
the states which predominantly contribute to transport. 
They can be read off from the master equation Eq. ^ 
and are summarized in Figure [5) 

Due to the strongly different coupling strengths, fill- 
ing the dot is much slower than emptying. Therefore one 
could expect that the double dot is predominantly empty, 
while single and double occupation are strongly sup- 
pressed. However, the isospin of an electron originating 
from the left lead has no overlap with the isospin polar- 
ization of the right lead. The direct transition |+)^ |0) 
is therefore forbidden and single occupation turns out to 
be more likely than an empty dot, as can be seen from 
the stationary occupation probabilities for a w — 1: 



/l-(a+l)\ 



4 

a + 1 




+ 0((a+l)2) (18) 



-f + ^ 1 +0((a + l)2). 



(19) 



with the isospin being parallel to n^. Due to the fact 
that filling the dot with a single electron occurs with rate 
4Tl and adding a second electron only with rate F^ the 



singly occupied state |+)^ is, in lowest order in 



1), 



four times more likely than an empty dot. Occupation 
of the triplet is even rarer than singlet occupation: it 
starts in order (a + 1)^, because it can only be reached 
via singlet occupation and subsequent decay to the state 

It is eminent from the flowchart Fig. [9] that there are 
several distinct cycles through which electrons are trans- 
ported from left to right: the transitions |+)^ ^ \S) and 
|0) ~^ \^) \'~)r ^^'^ then back to |0), or sev- 

eral sub-cycles via |T). As these cycles transfer electrons 
at different mean currents and with different statistics, 
it is clear that a complicated telegraph effect will lead to 
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increased noise. In contrast to the channel exclusion de- 
scribed in the previous sections this effect is not related 
to separated Hilbert spaces. 

V. IMPORTANCE OF SPIN 

We would like to remark that the divergence of the 
normalized cumulants for ~ 2T:m depends crucially on 
the inclusion of spin in the description of the system. 
The literature knows a number of examples where finite 
noise was found in similar, although not equal, double 
dot systems with spinless electrons There is also 

one example where infinite noise is predicted, although 
spin was not included.— At the end of this section we 



will discuss the relation to the present model. 

Neglecting spin in our model reduces the dimension- 
ality of the Hilbert space to four (instead of nine) since 
the distinction between singlet and triplet becomes im- 
possible. Instead there is only one doubly-occupied state, 
which is sufficiently described by its occupation probabil- 
ity p2- There are thus in principle 16 independent density 
matrix elements (instead of 81) of which only 6 (instead 
of 7) are independent. Again, we arrange the density ma- 
trix elements in a vector tt — (po,Pi,P2, 1), so that the 
master equation can be written in matrix form. Trans- 
forming to a new basis Sir — {po,pi/2 + l- nL,P2,Pi/2 — 
I • n^,! • (n^ X nL),Iz) similar to the above, the kernel 
of the spinless system W^i again assumes block-diagonal 
form for (f) = 27rm: 



-2ri 


Gl 
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e**2rflSin| 







-Gj^ - Gj^ 


Gr 





rsin| 








Gl 




Gl 
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Gl 
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-e'iVL sin| 


-irsini 


-e'TTflSinl 


-irsini 


-Tl - Tr 




















-Tl - Tr , 



(20) 



where, as above, — e*2 r,.(l ± cos ^) was defined. 

The block structure can be understood by realizing 
that, for an AB-phase = 27rm, charging the empty 
double dot |0) from the source always results in the sym- 
metric state |-f) = (|1,0) + |0, l))/\/2. From this state 
the electron may leave to the drain, resulting again in |0). 
In contrast to the spinful case, the symmetric state |+) 
cannot be charged with a second electron, so that |0) and 
|-|-) constitute a decoupled set of states, whose motion is 
described by the upper left block of Eq. (|20p . On the 
other hand, the doubly occupied state |1,1) may loose 
one electron to the drain, resulting in the antisymmetric 
combination |— ) = (|1, 0) — |0, 1))/V2, which can also be 
charged again from the source, but cannot be discharged 
to the drain. The two states {|— ) , |1, 1)} therefore also 
form a decoupled set, which is described by the mid- 
dle block of Eq. (pUj) . The remaining components of the 
isospin are unoccupied. 

In contrast to the situation with spinful electrons the 
statistics of the two subspaces arc the same, regardless of 
the coupling strengths. This is owed to the fact that both 
subspaces are two-dimensional and describe a single level, 
the statistics of which is symmetric in source and drain 
coupling. Correspondingly the statistics of the spinless 
model becomes that of two independent, non- interacting 
levels^i^ for = 27rm (see Fig. [TO]). 

The figure also shows that at = (2m -f l)7r the 
statistics becomes Poissonian. This is due to an isospin- 
blockade: Adding one electron to the empty system re- 



sults in the state (|1, 0) — « |0, 1))/a/2. This state cannot 
decay to the drain. On the other hand, the doubly occu- 
pied state may loose one electron, resulting in the com- 
bination (|1,0) + 1 10, l))/-\/2, which in turn cannot be 
refilled from the source. As a consequence, the system is 
trapped in the singly occupied state and transport events 
become increasingly rare as the flux approaches 27rm. re- 
sulting in Poissonian statistics. 

In summary, no super-Poissonian noise is predicted for 
any value of the magnetic flux when spin is neglected. In 




AB -Phase ^ 

FIG. 10: The normalized nth cumulants Knj n\ show no diver- 
gence for (\> = 27rm if spin is neglected. Instead they assume 
the values expected for a non-interacting two-level system. At 
= (2m -1- l)7r the statistics becomes Poissonian due to an 
isospin blockade effect. 
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Ref. diverging noise was also reported for spinless elec- 
trons as the flux approaches multiples of the flux quan- 
tum, (f) ~ 27rm. There a double dot Aharonov-Bohm 
interferometer with non-degenerate energy levels (level 
splitting Ae) was considered, with strong intra-dot in- 
teraction, so that double occupation is forbidden and a 
charge blockade mechanism can interrupt transport. It 
was pointed out that it is then important to perform the 
limit correctly as the system approaches the degeneracy 
point, i.e., taking first </) ^ and then Ae 0. Since we 
do not assume Coulomb interaction between the dots a 
charge blockade mechanism is not effective and the order 
of limits is uncritical. 

The spinless double-dot system with degenerate lev- 
els discussed here can be mapped onto a non-interacting 
quantum-dot spin-valvo^, with perfect lead polarization. 
The statistics of quantum-dot spin-valves have been an- 
alyzed in more general contexts, both without and with 
Coulomb interaction.— 



VI. CONCLUSION 

We have analyzed the full counting statistics of elec- 
tronic transport through an Aharonov-Bohm interferom- 
eter with two quantum dots embedded in its arms in 
the shot-noise regime. We found that for values of the 
Aharonov-Bohm flux that are integer multiples of the 
flux quantum, the second and higher cumulants diverge. 
This divergence is related to a separation of the Hilbert 
space of the double dot into disconnected subspaces that 
contain the spin singlet and triplet state for double occu- 
pancy, respectively. As the two subspaces have different 
transport statistics, the system exhibits electron bunch- 
ing that results in strongly super-Poissonian statistics. 
The coupling between the two subspaces, and therefore 
also the appearance of the divergence, is tunable by the 
Aharonov-Bohm flux. The inclusion of spin in the de- 
scription of the system has been shown to be crucial for 
the occurrence of diverging moments. Furthermore, we 
discussed how the divergence is cut off by the influence 
of an environment that relaxes either the isospin or the 
spin coherence. 
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APPENDIX A: PERTURBATIVE EXPANSION 
FOR LOWER ORDER CUMULANTS 

As mentioned in Sec. IIVI we make use of a recursive 
scheme to calculate the cumulants based on a master 
equation, which was developed by Flindt et al.r^ It is 
based on the assumption that the system's behavior is 
dominated by a single pole of the kernel W(x) = W(0)-f 
W(x). The cumulant generating function is then found 
to be 

Six) - (0| W(x) [1 + R(x)W(x)] |0) , (Al) 

where (0| and |0) are the left and right nuUvec- 
tors of W(0) and R(x) is the pseudoinverse R(x) = 
Q [W(x) - Six)]-' Q, with Q = 1 - |0) (0|. 

The cumulants can be obtained from Eq. (|A1[) using 
Brillouin-Wigner perturbation theory in x- This requires 
knowledge of the derivatives of the kernel. The first two 
cumulants (current and noise) are: 

Ki ^ i(0|w'|0) (A2) 

i 

K2 = 4(0|w"-2W'rw'|0), (A3) 

where the prime denotes a derivative w.r.t. the counting 
field X- 

The multiplication with left and right eigenvectors 
automatically picks out the correct eigenvalue. In 
our system we find the left nuUvector to be (0| = 
(1, 1, 1, 1, 0, 0, 0), so that it also automatically takes care 
of getting rid of the unneeded off-diagonal elements men- 
tioned above. 
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